Electronic energy level alignment at metal-molecule interfaces 

with a GW approach 
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Using density functional theory and many-body perturbation theory within a GW approximation, 
we calculate the electronic structure of a metal-molecule interface consisting of benzene diamine 
(BDA) adsorbed on Au(lll). Through direct comparison with photoemission data, we show that 
a conventional Go Wo approach can underestimate the energy of the adsorbed molecular resonance 
relative to the Au Fermi level by up to 0.8 eV. The source of this discrepancy is twofold: a 0.7 
eV underestimate of the gas phase ionization energy (IE), and a 0.2 eV overestimate of the Au 
work function. Refinements to self-energy calculations within the GW framework that account for 
deviations in both the Au work function and BDA gas-phase IE can result in an interfacial electronic 
level alignment in quantitative agreement with experiment. 
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There is considerable interest in using organic mate- 
rials as components in nanoscale energy conversion ap- 
plications, and thus a critical need has emerged for im- 
proved knowledge and control of the electronic struc- 
ture of metal-molecule interfaces. In particular, under- 
standing how molecular addition/removal energies (ion- 
ization energy, IE; and electron affinity, EA) are al- 
tered at a metal contact is fundamental to molecular- 
scale transport [IH3], energy conversion in organic pho- 
tovoltaics [4j [5], and photo- and electrocatalytic sys- 
tems [6]. 

Understanding metal-molecule interface electronic 
structure with spectroscopic accuracy poses significant 
challenges to standard first-principles approaches. Im- 
portant physical factors influencing electronic level align- 
ment include the magnitude of the interface dipole 
formed upon adsorption, molecular level broadening via 
hybridization with substrate states, and surface polar- 
ization effects on electron addition and removal energet- 
ics. While density functional theory (DFT) approaches 
within standard local and semi-local approximations can 
often describe interface dipoles [7HTU], hybridization, and 
work functions with good accuracy, prior studies [TT1 - 
[TT] have established that the impact of substrate po- 
larization, a non-local correlation effect, is absent from 
mean-field Kohn-Sham states. Self-energy corrections 
within the GW approximation can capture this effect, 
with a significant impact on gaps of adsorbed molecules 
(>1 eV for small aromatic molecules). GW methods 
can also significantly improve the IE and EA of gas- 
phase molecules [HI [18] compared to canonical semi-local 
Kohn-Sham DFT, where the energy difference between 
highest occupied molecular orbital (HOMO) and lowest 
unoccupied molecular orbital (LUMO) energies is under- 
estimated relative to the fundamental gap (i.e. IE - EA), 
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FIG. 1: Electronic density of states (DOS) of 1,4- 
benzenediamine (BDA) above the Au(lll) atop site. The 
inset shows geometry of the 4x4 supercell in color and periodic 
images in gray. The position of the highest occupied molec- 
ular orbital within DFT is -0.25 eV below the Fermi level. 
GoWo (this work) differs from the experimental measurement 
by 0.8 eV. The source of this discrepancy is twofold: a 0.7 eV 
underestimate of the gas phase ionization energy (IE), and 
a 0.2 eV overestimate of the Au work function. Accounting 
for these errors with a post-hoc correction, i.e. a rigid shift, 
produces a value (GW*) in excellent agreement with exper- 
iment. Results from resonant and ultraviolet photoemission 
spectroscopy (RESPES/UPS) [7] and their uncertainties are 
indicated by the shaded regions. 



even for the hypothetical "exact" exchange-correlation 
potential [T9U2T] . DFT frontier orbital energy differences 
can, however, provide accurate fundamental gaps if a ju- 
dicious approximation within a generalized Kohn-Sham 
framework is used [22] . 

Previous GW calculations of energy level alignment at 
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interfaces [23] [24j suggest significant improvement over 
DFT-GGA. In this work, we calculate the energy level 
alignment (shomo - £ Fermi, see Fig. [I]) at a proto- 
type metal-molecule interface, benzenediamine (BDA) on 
Au(lll), comparing our GW calculations directly with 
photoemission spectroscopy (PES) measurements. We 
find that modest inaccuracies in the constituent gas- 
phase BDA IE and Au(lll) work function within a stan- 
dard GoWo approach - using a plane-wave basis set and 
plasmon-pole models, and requiring sums over unoccu- 
pied single-particle states - are additive for the adsor- 
bate system in this case, leading to a discrepancy with 
measurements for the HOMO resonance of up to 0.8 eV. 
Refinements of self-energies within the GW framework 
that ameliorate deviations for the isolated constituents 
can lead to a predicted HOMO resonance in agreement 
with photoemission spectroscopy. 

GW calculations of metal-molecule interfaces pose sev- 
eral computational challenges. First, accurate evaluation 
of the Fock exchange requires explicit treatment of the 
semicore electrons, imposing the correct nodal structure 
on d states but leading to higher cutoffs and a need to 
treat more electrons [25 j. Second, metals often require 
a dense k-point sampling and the relevance of plasmon 
pole approximations can be questionable [25j[26]. Third, 
for a hybrid interface comprised of a molecule (with lo- 
calized states) and a metal (with delocalized states), self- 
energies for the constituent systems are of very different 
magnitudes, and taking the Kohn-Sham eigenstates as 
the quasiparticle wavefunctions may no longer be a good 
approximation [27]. Fourth, for level alignment between 
states of disparate character, we require absolute conver- 
gence, and, in this case, since Au self-energy corrections 
converge differently than those for BDA states, a large 
number of unoccupied states (N c ), together with a good 
extrapolation scheme [28] [29] , can be necessary. This is 
a particular challenge given the concomitant need for a 
large supercell. 

Our GW calculations are performed using the Berke- 
leyGW [30] and Abinit [31] codes, following a well- 
established Go Wo approach [32 . Equilibrium geome- 
tries of molecular BDA in the gas-phase and physisorbed 
on Au(lll) are determined using DFT [33] within the 
PBE [34 generalized gradient approximation (GGA). 
The molecule is flat relative to the surface, at a height of 
3.5 A above the topmost layer of Au. This is consistent 
with relaxed geometries obtained using a van der Waals 
corrected density functional [35] . Norm- conserving pseu- 
dopotentials [20] are used with a plane- wave basis (60 
Ry cutoff) for structural relaxations and includes 5s and 
bp semi-core states for Au. The surface is modeled with a 
4x4 supercell containing 4 layers of Au (roughly 9 x k^ rrni , 
where /cpermi is the Thomas-Fermi wavevector of gold), a 
single BDA molecule (see inset Figure [TJ, and the equiv- 
alent of 10 layers of vacuum. The metal work function 
and magnitude of I^homo — £ Fermi | change less than 0.05 



eV in DFT when the depth of the slab is doubled. The 
theoretical in-plane bulk lattice parameter is used for Au 
(a=4.18 A). The supercell Brillouin zone (BZ) is sampled 
using a 4x4x1 k-grid. Gas-phase BDA is modeled using 
the same supercell, in the absence of Au, and using a 
Coulomb truncation. 

For all calculations, 6 Ry planewave expansion cutoff 
is used for the dielectric function, which, for the major- 
ity of our work, is extended to finite frequencies with a 
generalized plasmon pole (GPP) model [32] . Doubling 
this cutoff and the number of unoccupied states used in 
constructing the response function results in negligible 
changes to BDA gas-phase IE, ~0.15 eV. Updates to G 
and W use quasiparticle energies from the previous cycle 
and a linear fit to a coarse sampling of self-energy cor- 
rections to high energy states (s > 6 eV above vacuum). 
The number of states used to construct e _1 was held 
fixed at 2048 bands. For BDA-Au(lll) interface calcu- 
lations, our sum over the unoccupied subspace includes 
more than 1400 conduction bands (30 eV above ^Fermi), 
a number which, as we show, still falls considerably short 
of convergence. Gas-phase results are based on a sum of 
over 5100 conduction bands (~ 80 eV above vacuum) for 
the Coulomb-Hole term. For calculations of bulk Au, we 
compare with results from the Godby-Needs [36 GPP 
model, as well as a an explicit evaluation of e _1 (q, uo) for 
more than 200 frequencies up to about 100 eV [37 . For 
bulk Au, the BZ is sampled with a 14 3 k-point grid, and 
500 conduction bands is used (~ 600 eV above ^Fermi)- 

Table [J summarizes our results for the IE and EA of 
gas-phase BDA at different levels of theory. Relative 
to measured photoemission, DFT-GGA underestimates 
the IE by over 3 eV, consistent with previous work pQ. 
Our GW calculations are a significant improvement over 
DFT-PBE, within 0.7 eV or better of experiment (after 
extrapolating the unoccupied states to infinity), depend- 
ing on whether G or W is updated. The error relative to 
experiment for GiWi is 0.5 eV, or just 7%. Remarkably, 
with just 1000 unoccupied states and no extrapolation, 
the IE is about 1 eV smaller, illustrating the slow conver- 
gence of the IE with respect to unoccupied states. Dou- 
bling the number of unoccupied states reduces the IE by 
only 0.2 eV. Use of an extrapolation scheme for this slow 
convergence, associated with the Coloumb-Hole term of 
the GW self-energy, is therefore crucial when comparing 
GW to experiment. Extrapolations 7V C — »• oo are de- 
termined by fitting the Coulomb-Hole term to the form 
i 

a + N c x * , where x* is determined from a similar fit to 
the static Coulomb hole screened-exchange (COHSEX) 
approximation (using the same set of convergence param- 
eters). Within static COHSEX, the asymptote a does not 
directly depend on 7V C , allowing for clean optimization of 
the exponent. Encouragingly, our scheme seems consis- 
tent with results computed using a completion method 
recently proposed by Deslippe et al. [29] . 
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TABLE I: Energy levels for the frontier molecular orbitals in 
the gas phase calculated using different electronic structure 
methods. All values are in eV. Results for DFT-PBE cor- 
respond to Kohn-Sham eigenvalues of the neutral molecule. 
Values for N c — >• oo are based on an extrapolation of the 
Coulomb- Hole term (see text). Uncertainties in the experi- 
mental ESI EH value of the vertical IE is ± 0.03 eV. 



For bulk Au, we find that the GPP models of 
Hybertsen-Louie and Godby-Needs leave the DFT-PBE 
bandwidth relatively unchanged, in adequate agreement 
with XPS measurements [40] (overestimate of 12% and 
9% respectively). Using a fully frequency dependent di- 
electric function results in a 6% overestimate relative to 
experiment. Better agreement will likely require going 
beyond the RPA gTJ[42]. Interestingly, DFT-PBE pro- 
vides the best Au(lll) work function [43], 5.2 eV (within 
0.05 eV of experiment [44 ). With either of the GPP val- 
ues used here, the Au(lll) work function is larger than 
experiment by 0.5 eV, and by 0.2 eV with a fully fre- 
quency dependent dielectric function. For comparison, 
Faleev et al. [45 report a Go Wo (and numerical e _1 ) 
work function for Al(lll) which is 0.06 eV too small com- 
pared to experiment (4.18 eV vs 4.24 ± 0.02 eV), with a 
self-consistent GW approach producing nearly identical 
values. 

To identify the HOMO energy of the BDA adsorbate, 
we project the DFT Kohn-Sham Hamiltonian, H K , of 
the BDA-Au(lll) supercell onto the DFT-PBE orbitals 
calculated from the isolated molecule, |^). For BDA ph- 
ysisorbed on Au(lll), |(^homo|^)| 2 ^ 0.9 at the T point, 
indicating such a projection is spectroscopically mean- 
ingful in this case. We note that as the BDA HOMO 
resonance is not the highest occupied state in our calcu- 
lation, in general we would not expect its value to agree 
with UPS, at least within the semi- local KS framework 
used here. Indeed, evaluating (lhomo\H k \lhomo) > we 
obtain -0.25 eV, a resonance value too shallow compared 
to photoemission, which place the HOMO -1.4 ± 0.1 eV 

below Spermi- 

To evaluate GW self-energy corrections for the BDA 
adsorbate HOMO energy, we follow the above approach 
and evaluate the matrix element (^homo|S k |^homo)- 
Since E K is approximately diagonal in the molecular ba- 
sis {|^j)}, this approach is approximately equivalent to 
evaluating diagonal and off-diagonal elements of the full 
energy-dependent self-energy matrix, diagonalizing, and 
projecting onto the surface states to identify the adsor- 
bate HOMO resonance energy measured spectroscopi- 
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FIG. 2: Scon for the isolated (dashed line) molecule and 
molecular resonance (solid line) in the monolayer as a func- 
tion of the number of conduction bands used in the evalua- 
tion of E. Although the absolute value of these terms con- 
verge slowly, by N c = 600 their difference, AE CO rr (inset), has 
stabilized to 1.7 eV, attributable to non-local static correla- 
tions from the metal substrate. The nonlinear behavior of 
S CO rr(A/ r c ) for N c < 200 reflects the character of low energy 
conduction bands, which are highly system specific. 



cally [T7JE2]- This workaround, with its substantially- 
reduced computational cost, is strictly valid in a "weak- 
coupling" limit, where both H K and E^ are diagonal in 
the basis of gas-phase orbitals, as is the case here. 

To understand the adsorbate result relative to the gas- 
phase, we follow previous work [TT] and partition the 
self-energy correction into two contributions: the Fock 
exchange, Ex; and the portion containing static and dy- 
namical correlation, E C orr- We find that Ex f° r the ad- 
sorbate HOMO differs that obtained for gas-phase BDA, 
Ex (adsorbate) — Ex (gas phase) = 0.4 eV, which can be 
understood in terms of the non-zero overlap of the BDA 
HOMO with the Au wavefunctions. Unlike Ex? S corr 
involves a difficult-to-converge sum over the unoccupied 
space (as with gas-phase BDA). However, from Fig.[2j the 
difference between E corr of the isolated molecule and the 
adsorbate monolayer converges much faster, with a mod- 
est sum of 600 bands (Fig. [2| to 1.7 eV. As with benzene 
on graphite [TT], this response is due almost entirely to 
static polarization: an electrostatic image charge model 
(with a calculated [46] image plane of 1.47 A above the 
Au surface) predicts a value of AE corr = 1.8 eV. 

As shown in Fig.[T] our Go Wo corrections, in the limit 
of N c — )> oo, lead to a BDA adsorbate HOMO energy of 
-0.64 eV, a significant underestimate of the experimental 
value of -1.4 eV. Given the converged self-energy correc- 
tions for the isolated molecule, the metal work function, 
and the change in the correlation energy upon adsorp- 
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tion, we can understand the disagreement between theory 
and experiment for the adsorbate system as originating 
with the underestimate of the gas-phase BDA IE and the 
overestimate of the Au(lll) work function. If we account 
for the discrepancies of the isolated systems with a rigid 
shift (GW* in Fig. [TJ, good agreement between theory 
and experiment for the composite system is obtained. 

Our results illustrate that the accuracy of energy level 
alignment at a metal-molecule interface with a given GW 
approach is limited by its ability to describe the IE of 
an isolated molecule and the metal work function. In a 
stronger-coupling limit, dynamical contributions to elec- 
trode polarization would become important, as has been 
noted before [TTJ [131 [14] . In such a case, classical static 
polarization models are less valid, and S C orr obtained 
from GW must be used. Furthermore, S K may no longer 
be diagonal in the basis of gas phase orbitals, necessitat- 
ing full evaluation of the self-energy operator. 

In conclusion, through direct comparison with pho- 
toemission spectroscopy, we have demonstrated the ad- 
vantages and limitations of an existing Go Wo approach 
in describing the electronic structure of a molecule ad- 
sorbed to a metal substrate. GW improves upon PBE, 
particularly in its inclusion of nonlocal correlation, but 
is limited, at least in the approach considered here, by 
its ability to predict gas-phase IEs and work functions. 
Computationally-tractable refinements that improve ac- 
curacy of these quantities will result in better quantita- 
tive agreement for level alignment. 
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